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The lightest supersymmetric particle may decay with branching ratios that correlate with neu- 
trino oscillation parameters. In this case the CERN Large Hadron Collider (LHC) has the potential 
to probe the atmospheric neutrino mixing angle with sensitivity competitive to its low-energy de- 
termination by underground experiments. Under realistic detection assumptions, we identify the 
necessary conditions for the experiments at CERN's LHC to probe the simplest scenario for neutrino 
masses induced by minimal supergravity with bilinear R parity violation. 

PACS numbers: 11.30.Pb,12.60.Jv,14.60.Pq,95.30.Cq 



I. INTRODUCTION 



The CERN Large Hadron Collider (LHC) will provide high enough center-of-mass energy to probe directly the 
weak scale and the origin of mass [lj-Q . In addition to its designed potential, here we show how LHC searches for new 
physics at the TcV region may provide an unexpected opportunity to probe neutrino properties, currently determined 
only in neutrino oscillation experiments Q, shedding light on some of the key issues in neutrino physics. We illustrate 
how this works in a class of supersymmetric models where the lepton number is broken, together with the so-called 
R parity symmetry Even when the latter holds as a symmetry at the Lagrangian level, as in some SO(10) 
unification schemes, R parity breaking may be driven spontaneously by a nonzero vacuum expectation value of an 



SU(3) ® SU(2) ® U(l) singlet sneutrino [9H12| . In this case the low-energy theory is no longer described by the minimal 



supersymmetric standard model, but contains new R parity violating interactions 
of this scenario leads to an effective model with bilinear violation of R parity 



134l5|. The simplest realization 
16h 2011 . The latter constitutes the 



minimal way to break R parity in the minimal su persym metric standard model and provides the simplest intrinsically 



supersymmetric way to induce neutrino masses [2ll424| . Its main feature is that it relates lightest supersymmetric 



particle (LSP) decay properties and neutrino mixing angles 



25-27 



Here we demonstrate that indeed, under realistic assumptions, the simplest scenario for neutrino masses in su- 
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persymmetry (SUSY) with bilinear violation of R parity can be tested at the LHC in a crucial way and potentially 
falsified. We identify the regions of minimal supergravity (mSUGRA) parameters, event reconstruction efficiencies and 
luminosities where the LHC will be able to probe the atmospheric neutrino mixing angle with sensitivity competitive 
to its low-energy determination by underground experiments, both for 7 and 14 TeV center-of-mass energies. 

For the sake of definiteness, we consider the minimal supergravity model supplemented with bilinear R parity 
breaking 22H24j added at the electroweak scale; we refer to this scenario as RmSUGRA. In this effective model one 



typically finds that the atmospheric scale is generated at tree level by a weak-scale neutralino-cxchange seesaw, while 
the solar scale is induced radiatively [3] . The LSP lacks a symmetry to render it stable a nd, given the neutrino mass 



scales indicated by oscillation experiments, typically decays inside the LHC detectors [22|, |23j, |25[ 1 . As an illustration 
we depict the neutralino LSP decay length in Fig. [TJ We can see from Fig. Qlthat the expected decay lengths are large 
enough to be experimentally resolved, leading to displaced vertex events [33l . 134 1 . 
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Figure 1: \\ decay length in the plane mo,mi/ 2 for Ao = —100 GeV, tan/3 = 10 and jj, > 0. 



More strikingly, one finds that in such a RmSUGRA model one has a strict correlation between neutralino de- 
cay properties measurable at high-energy collider experiments and neutrino mixing angles determined in low-energy 
neutrino oscillation experiments, that is 



BR{x1 -> m ± TU =f ) 
BR{xl -> rilUT) 



tan" atm ~ — £i ^ . (1) 



The derivation of Eq. ([T]) can be found in 25|. In short, the relation between the neutralino decay branching ratio 
and the low-energy neutrino angle in the bilinear model can be understood in the following way. At tree-level in 
RmSUGRA the neutrino mass matrix is given by [22I 




eff = a A^±l A A \ A eA, A, A J (2) 



where Aj = fiVi+vo^i and and i>j are the bilinear superpotential parameters and scalar neutrino vacuum expectation 
value, respectively. Equation ([2} is diagonalized by two angles; the relevant one for this discussion is the angle 
tan 023 = — -j^-- One can understand this tree-level mass as a seesaw- type neutrino mass with the right-handed 
neutrino and the Yukawa couplings of the ordinary seesaw replaced by the neutralinos of the minimal supersymmetric 



1 We may add, parenthetically, that such schemes require a different type of dark matter particle, such as the axion [28l |. Variants with 
other forms of supersymmetric dark matter, such as the gravitino I29H32I1 . are also possible. 
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standard model and couplings of the form cAj, where c is some combination of (generation independent) parameters. 
These couplings, which determine (the generation structure of) the neutrino mass matrix, also determine the couplings 



If - WT and xl 



25j . Taking the ratio of decays to different generations the prefactors c drop out 
and one finds Eq. (TTJ), when the angle tan #23 is identified with the atmospheric neutrino angle. One- loop corrections 
tend to modify this relation, but, as long as the loop corrections are smaller than the tree-level neutrino mass, Eq. (TTJ) 
is a good approximation [251 ] . 

In other words, as seen in Fig. [5J the LSP decay pattern is predicted by the low-energy measurement of the 
atmospheric angle 21, 25 1. currently determined by underground low-energy neutrino experiments Q, 
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Figure 2: Ratio of y± decay branching ratios, Br(x? — ► ^q'q) over Br(x? — > Tq'q) in terms of the atmospheric angle in bilinear 
R parity violation [25| . The shaded bands include the variation of the model parameters in such a way that the neutrino masses 
and mixing angles fit the required values within 3a. 



In this paper we show how a high-energy measurement of LSP decay branching ratios at the LHC allows for a 
redetermination of # a tm and hence a clear test of the model. We provide quantitative estimates of how well this ratio 
of branchings should be measured at LHC in order to be competitive with current oscillation measurements. This 
issue has already been addressed but only at the parton level, using some semirealistic acceptance and reconstruction 
cuts, and for just one specific mSUGRA point [35 1. 



II. FRAMEWORK OF OUR ANALYSIS 



Our goal is to present a more detailed analysis of the LHC potential to measure the LSP branching ratios required 
to test the relation shown in Eq. (fTJ , going beyond the approximations made in the previous work of Rcf. [35j . The 
generation of the sup ersymmetric spectrum and decays in the scope of the RmSUGRA model was carried out using the 
SPheno package 36j 2 . The event generation was done employing PYTHIA 13 71 with the RmSUGRA particle properties 
being passed into it in the SUSY Les Houches accord (SLHA) format 38|, |39|. Jets were defined using the subroutine 
PYCELL with a cone size of AR = 0.4. 



2 An updated version including bilinear R parity violation can be obtained at http : //www . physik . uni-wuerzburg . de/~porod/SPheno . html. 
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A striking property of RmSUGRA models is the existence of displaced vertices associated to the LSP decay [34 1 . 
We use the detached vertices to probe the LSP branching ratio relation Eq. ([TJ . In order to mimic the LHC potential 
to study displaced vertices we use a toy detector based on the ATLAS technical proposal Q . 

We begin our analysis demanding that the events pass some basic requirements to guarantee that they will be 
triggered by the experimental collaborations. This is done because the LHC experiments have not defined so far 
any specific strategy to trigger displaced vertices with such high invariant mass, therefore, we restricted our analysis 
to events that would be accepted by the ongoing analyses. We accept events passing at least one of the following 
requirements, denoted as cut CI, 

1. the event has one isolated electron or a photon with pt > 20 GeV; 

2. the event has one isolated muon with px > 6 GeV; 

3. the event has two isolated electrons or photons with px > 15 GeV; 

4. the event has one jet with pr > 100 GeV; 

5. the event has missing transversal energy in excess of 100 GeV. 

Next, in cut C2, we require that at least one of the neutralinos in the event decays beyond the primary vertex 



point, that is, outside an ellipsoid 34 1 



y \ , i z 



1, (3) 



where the z axis is taken along the beam direction. We made a conservative assumption, since we are not performing 
a detailed detector simulation, that the ellipsoid dimensions are 5 times the ATLAS expected resolutions in the 
transverse plane (6 xy = 20 fiva) and in the beam direction (S z = 500 /an), in order to ensure that the neutralino 
displaced vertex is distant of the primary vertex. We also demand that all tracks must be initiated inside the pixel 
inner detector within a radius of 550 mm and z axis length of 800 mm. A detached vertex complying with these 
requirements we called signal vertex. 

In order to check relation Eq. (fTJ we looked for detached vertices presenting a W associated to them and we must 
isolate the LSP decays into W/i and Wt. Moreover we consider only hadronic final states of the W as a necessary 
condition for the identification of the lepton flavor. In cut C3, which is designed for the W reconstruction, we 
require two jets with charged tracks intersecting the neutralino resolution ellipsoid, and invariant mass between 60 
and 100 GeV. In order to be sure that the W reconstruction is clean, we further impose that the axes of other jets 
of the event to be outside of a cone AR = 0.8 of the W jets' axes. Note that this cut should eliminate standard 
model (SM) backgrounds coming from displaced vertices associated to 6's or r's. To guarantee a high quality in the 
reconstruction of the displaced vertices we impose that the W decay jets must be central, having pseudorapidities 
\r)\ < 2.5; this constitutes our cut C4. The events passing the above requirements most probably originate from LSP 
decay, having basically no sizable standard model background, except for instrumental backgrounds and beam-gas 
interactions. 

A signal vertex is classified as originating from the LSP decay into a /iW pair if it presents a fi^ and a hadronically 
decaying W stemming from the displaced vertex with transverse momentum px > 6 GeV and I77I < 2.5. In the r 
case we demanded that the t ± associated to a detached W possesses pr > 20 GeV and \r]\ < 2.5. These requirements 
are called C5. 

Detecting taus is somewhat more complicated than detecting muons, so one needs to be more careful in recon- 
structing the tW pair displaced vertex. The following criteria, denoted C6, arc used to separate the detached vertices 
exhibiting a through its 1- and 3-prong decay modes. We check also that the secondary displaced vertex from 
tau decay does not spoil the signal vertex; i.e., we verify that the tau decay products point towards the LSP decay 
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vertex within the experimental resolution. We define the neutralino resolution ellipsoid as the ellipsoid centered at 
the displaced vertex position of neutralino, Vx, with axes 5 xy = 12 /mi and 8 Z = 77 fim based on ref. [3J. Let p pr0 ng be 
the momentum of either 1-prong tau decay or the sum of momenta of the 3-prong decays. Let also v 2 be the position 
of the secondary vertex coming from r. We verify whether the line along p pr0 ng, crossing v 2 intersects the neutralino 
resolution ellipsoid. For this we require that for each t, the discriminant of quadratic equation for parameter t 



Pprong^ + V 2 




(4) 



be equal to or greater than zero. In previous 35[ analysis only 3-prong tau decays modes were considered. 



An additional cut C7 was applied to 3-prong tau events, i.e. we also require that one of the prongs has a transverse 
momentum pt > 9 GeV while the other two have px > 2 GeV. In addition we check if all prongs lie within a cone 
radius of AR < 0.2 around the tau direction obtained from the prongs' tracks. 

Finally we require that the signal lepton (/x or r) be isolated; cut C8. \x isolation demands that there are no 
other tracks whose total transverse energy satisfies Et > 5 GeV within a cone AR > 0.3. The t was required to 
be isolated using the same criteria as for the muon, but for an annulus of outer radius AR — 0.4 and inner radius 
AR = 0.1. Isolation of the leptons is a needed requirement to eliminate events presenting lcptons generated inside 
jets and constitutes an important cut to reduce potential backgrounds. 



III. RESULTS AND DISCUSSION 



In order to access the effects of the above defined cuts C1-C8 we present detailed information on their effects 
for the mSUGRA SPSla benchmark point characterized by m 1/2 = 250 GeV, m = 100 GeV, A = -100 GeV, 



tan /? = 10, and sgn(/i) = +1. This allows us to compare our results with the one previously obtained in [351 ] . For the 



default solution of SPheno to the neutrino masses and mixings, the relevant neutralino branching ratios are 

BR(x° -> W ± n T ) = 5.4% BR(x? -> W ± t^) = 6.2% BR(x? -> Zv) = 1.2% 

BR(x? e ± r^v) = 11.5% BR{x" -> fj^i^v) = 24.3% BR(x? -> r ± r T t/) = 36.4% 

BR(x? — » bbv) = 14.7%; (5) 

with the R parity parameters being 

ex = 0.0405 GeV, e 2 = -0.0590 GeV, e 3 = 0.0506 GeV, 

vx = -0.0027 GeV, v 2 = 0.0042 GeV, v 3 = -0.0033 GeV. 

Furthermore, for this choice of parameters the neutralino decay length is cr = 1.1mm, and it travels an average of 
4.4 mm in the laboratory. 

From Table U we see that the vast majority of the events pass the trigger requirements CI, as expected. For the 
SPSla SUSY point, the LSP decay length is sufficiently long to guarantee that a sizeable fraction of its decays take 
place away from the primary vertex; this reflects as a high efficiency for passing the cut C2. We have focused our 
attention to events presenting a W decaying into two jets through C3. It is interesting to notice that 63% of the W 
hadronic decays are in the form of two jets. Additional suppression of the signal by C3 comes from the matching of 
the sum of momenta of the charged tracks pointing to the detached vertex and the jets reconstructed using PYTHIA. 
To further illustrate the W decay, we present in Fig. [3] the jet-jet invariant mass distribution. As we can see, this 
distribution is clearly peaked around the W mass and a good fraction of the two jets reconstructed as associated to 
the LSP decay pass the cut C3. The observed high efficiency of cut C4 shows that the W's produced in the LSP 
decay are rather central. 
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Table I: Fraction of events passing the successive cuts CI— C8 used for the event reconstruction at the SPSla mSUGRA point. 




Figure 3: From top to bottom: xl ~* jjX without cuts, with cut on lepton isolation (fi or r) and with all other cuts leaving 
free the invariant mass range. 

We also learn from Table Q] that detached vertices presenting a W possess around 60% of the time an energetic //* 
or t ± complying with C5. Moreover the cuts C6 and C7, which ensure the quality of the r reconstruction, reduce 
significantly the number of T4 A± r =F events. Finally the isolation cut C8 turns out to be quite important significantly 
reducing the signal. 

For the parameter point SPSla, the expected efficiencies for the reconstruction of [iW and tW decays are 0.107 
and 0.0098 respectively, where in the last we have added 1- and 3-prong hadronic decays. When the r decays into a /i 
and neutrinos, the event was computed as being a fiW decay if the [i passes the cuts. This was included appropriately 
in our calculations. Taking into account the total SUSY production cross section (41 pb) at 14 TeV, an integrated 
luminosity of 100 fb" 1 and these efficiencies we anticipate that the number of observed fiW and tW events after cuts 
to be 

Nf, =32000 jyhadron = 3333 

where N l ™ dlon = N^^iln + iV^ prong . Therefore, the statistical accuracy of the ratio R = BR{x x 
H ± W T )/BR(xi -> t ± VF =f ) is expected to be a(R)/R = Jl/N^ + 1/N T sa 0.015. In the case one takes into ac- 
count only the three-prong decays of the tau, as in Ref. [35[, the statistical error of this ratio increases to w 0.053. 
Moreover, as expected, there is a degradation of the accuracy in the determination of this ratio of branching ratios in 



a more realistic analysis; the result obtained in 35j is ~ 0.028 
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In the evaluation of the above efficiencies we have not taken into account multiple interactions at the LHC as needed 
for the high luminosity run. Therefore, we reevaluated the detection efficiencies for muons and taus with multiple 
interactions switched on in PYTHIA. We found that these efficiencies were only slightly degraded by the occurrence of 
pilcup, that is, we obtained that the efficiencies for muon reconstruction are reduced to 0.102 and for tau are 0.008 68 
in hadronic mode and 0.000 94 in the 3-prong mode. In our analyses we took into account the effect of multiple 
interactions. 

For the sake of comparison, we present a detailed analysis for a different mSUGRA point that is TO1/2 = 500 GeV, 
mo = 500 GeV, = — lOOGcV, tan/3 = 10, and sgn(/i) = +1. Once again using SPheno, we obtain that the 
neutralino branching ratios larger than 1% are: 

BR(x° -> W ± fi T ) = 22.9%, BR(x? -> W ± t t ) = 25.2%, BR(x? Zv) = 25.1%, 

BR(x? H- vh°) = 16.9%, BR(x? -4 r ± T =F i/) = 3.4%, BR(x? -> bbv) = 2.9%; 

and the corresponding R parity parameters are 

ei = 0.1507 GeV, e 2 = -0.1507 GeV, e 3 = 0.1507 GeV, 

v-l = -0.0056 GeV, v 2 = 0.0058 GeV, v 3 = -0.0054 GeV. 

As we can see, the neutralino LSP decays are dominated by the two-body ones, in contrast with the SPSla point 
where the three-body decays mediated by light scalars are dominant. Because of its heavier spectrum, the total 
SUSY production for this parameter point is smaller than the SPSla one; however, the cross section loss is partially 
compensated by the higher branching ratios into fiW and tW. 

The total cross section for this case is 832.0 fb and our analyses indicate that the reconstruction efficiency for fj,W 
decays is 0.203 while the tW decays are reconstructed with an efficiency of 0.035, where we did not take into account 
pilcup. The inclusion of this effect leads to a tiny reduction of the reconstruction efficiencies that become 0.199 for 
/iW and 0.033 for tW. On the other hand the efficiency for reconstructing a tW event in the 3-prong mode is 0.012. 
Notice that these efficiencies are larger for this mSUGRA point than for the SPSla because the neutralino is heavier 
and, consequently, its decay products are more energetic and pass the cuts more easily. The expected total number 
of reconstructed events after cuts for this SUSY point is 7V M = 5171 and 7V^ adron = 933 where we have included the 
pilcup effects. Therefore, the expected statistical error on the ratio R becomes ss 0.036, or « 0.056 when we only use 
3-prong taus as in [35)]. As we can sec, the statistical error on the ratio R increases as m 1( / 2 (LSP mass) increases due 
to the reduction of the SUSY production cross section despite the increase in the detection efficiencies. 

We evaluated the reconstruction efficiencies as a function of mo <£> m\n for Aq = —100 GeV, tan/3 = 10 and 
sgn(/i) = +1 and our results are depicted in Fig. S) As we can see from the left panel of this figure, the yW decays 
exhibit a high reconstruction efficiency, i.e., between 10% and 20%, in a large area of the parameter space, degrading 
only at large m^. On the other hand, the tW reconstruction (see right panel of Fig. [4]) is at most 3.5%, indicating 
that the statistical error on the ratio R is going to be dominated by these events. 

We present in Fig. [5] the attainable precision a(R)/R with which the correlation R can be measured as a function 
of mo <£> mi/2 for Aq = —100 GeV, tan/3 = 10, and sgn(^) = +1 for an integrated luminosity of 100 fb -1 and a 
center-of-mass energy of 14 TeV. We require in all plots that at least 5 events of reconstructed taus are observed. In 
the left panel of this figure we present the expected statistical error on the ratio R assuming no systematic errors on 
the determination of the reconstruction efficiencies, while in the right panel we consider a more conservative scenario, 
where we anticipate a systematic error of 10% in each of the reconstruction efficiencies. One can see from this panel 
that the precision drops as m^ 2 grows since the neutralino production rates from squark/gluino cascade decays also 
decrease with increasing m 1 / 2 values. Therefore, if the systematic errors of the efficiency determination are negligible 
the LHC collaborations should be able to probe with a very good precision (< 10%) the ratio R for m\/2 < 650 GeV, 
which correspond to an LSP mass up to ~ 270 GeV. The inclusion of systematic errors at the level assumed in the 
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Figure 4: Reconstruction efficiencies of fiW (left panel) and tW events (right panel) as a function of mo 8^1/2 for A) = —100 
GeV, tan/? = 10 and sgn(/x) = +1 including the effect of pileup. The red (dark shaded) area corresponds to the region where 
stau is the LSP, while the yellow (light shaded) area represents the region excluded by LEP. 

right panel of Fig. [5] increases the uncertainty in R; however, it is still possible to perform an accurate test of the 
RmSUGRA scenario. 

Note that in Fig. [5] we also present results for the 7 TeV run of the LHC. For this case one can see that the LHC has 
a much more limited capability of probing the ratio R, since the reach of this run covers only up to m±/2 J$ 300 GeV. 
Still, although large, the statistical errors in this region [0.3 < a{R)/R < 0.5], due mainly to the small anticipated 
integrated luminosity, which we have taken to be 1 fb -1 , allow a determination of the atmospheric angle comparable 
to that obtained at low energies. 

In the left panel of Fig. [6] we show the dependence of the attained precision as a function of the neutralino mass for 
luminosities of 2, 10, and 100 fb -1 . For small neutralino masses the SUSY production cross section is large enough to 
guarantee that the statistical errors are small; therefore, the uncertainty on the ratio R is dominated by the assumed 
systematic errors on the reconstruction efficiencies, even for an integrated luminosity of 2 fb -1 . As the accumulated 
luminosity increases the LHC experiments will be able to probe higher neutralino masses; however, the precision 
worsens due to the increase of statistical errors. We can also see clearly that increasing the luminosity allows a more 
precise measurement of R as expected. Moreover, one can probe LSP masses up to 250 (320 or 370) GeV for an 
integrated luminosity of 2 (10 or 100) fb -1 . 

From the right panel of Fig. [5] we estimate the luminosity needed to measure R with a given precision for several 
LSP masses. For instance, let us consider m^o = 250 GeV. In this case R can only be measured with a precision 
a(R)/R ~ 50% with 2 fb -1 , while this error can be brought down to 20%, i.e., close to the limit set by the systematic 
uncertainties, with 50 fb -1 . 



IV. CONCLUSIONS 

We have demonstrated how the LHC may have the potential of probing neutrino mixing angles with sensitivity 
competitive to their low-energy determination by oscillation experiments. This analysis was carried out, for the sake 
of concreteness, in the simplest scenario for neutrino masses induced by minimal supcrgravity with R parity violation 
as framework. In this class of models, the smoking gun for the neutrino mass generation mechanism is the ratio of 
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Figure 5: Precision in the determination of the ratio R in the plane myj x mo for a luminosity of 100 fb - , center-of-mass 
energy of 14 TeV, Ao = —100 GeV, tan/3 = 10, and sgn(^i) = +1. In the right (left) panel we did (not) include a possible 
systematic uncertainty in the extraction of the efficiencies for the channels /j,W and tW . The stars in the right panel represent 
the results for the 7 TeV run with an integrated luminosity of 1 fb _1 . The shaded areas represent the same as in Fig. [4] 
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Figure 6: The left panel displays the achievable precision in the ratio R as a function of the neutralino mass m^o for luminosities 
of 2, 10, and 100 fb -1 at 14 TeV whereas the right panel contains the foreseen statistical error on R as a function of the integrated 
luminosity for several LSP masses. 
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branching fractions of ncutralino decaying into fj,W and tW, as this fraction is related to the atmospheric neutrino 
mixing angle in RmSUGRA models. 

Under realistic detection assumptions we have made a detailed analysis of the reconstruction of ncutralino decays, 
as well as of the cuts needed to characterize the signal. After that we determined the attainable precision on the 
measurements of the ratio R given in Eq. (QJ. Comparing with a previous parton level study, we improved the 
reconstruction efficiencies of muons as well as taus. 

We showed that the 7 TcV run of the LHC will have a somewhat weak potential for probing the RmSUGRA 
model, since it is statistics limited. Still, precisions comparable to the low-energy determination should be reached. 
In contrast, a 14 TeV run with 100 fb _1 integrated luminosity will be able to probe a large fraction of the parameter 
space with a good precision, as seen in Fig. [5] In fact, our analyses suggest that the error on R will be dominated 
by the systematic ones on the reconstruction efficiencies of the decay \iW and tW, with the statistical errors being 
under control. 

In short, we find that in this case the atmospheric mixing angle may be probed relatively neatly. In fact, a 
determination of R within a given error translates into a prediction of the atmospheric mixing angle with an error 
of very similar size. Needless to say, what we have presented is only one example of a class of LSPs. There are 
other variant schemes based on alternative supersymmetry and/or R parity breaking, where other states emerge 



as LSP and similar correlations to other neutrino mixing angles appear [4lM43| . These would, however, require 
separate dedicated studies. We encourage the particle detector groups ATLAS and CMS to add the test of such 
possibilities to their physics agenda, as this might lead to a tantalizing synergy between high-energy accelerator and 
low-energy nonaccclcrator searches for new physics. Studies with the real LHC data may also make it possible to 
probe, at some level, the mass scale characterizing atmospheric neutrino oscillations, as well as the angle characterizing 
solar neutrino oscillations, an issue to be taken up separately. 
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